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1 Introduction 



As it is known, when dealing with the estimation of the integral Ia{f) — 
f{x)da{x), a{x) being a positive measure on [a,b] by means of an n-point 
Gauss-Christoffel quadrature rule, In{f) — J2]=iAjf{xj) such that Ia{P) — 
In{P) for any polynomial of degree up to 2n — 1, the effective computation 
of the nodes {xjY-^i and weights {Aj}"^^ in /n(/) has become an interesting 
matter of study both numerical and theoretical. As shown by Gautschi (see 
[16], [17] or [18]) among others, here the basic fact is the three-term recurrence 
relation satisfied by the sequence of orthogonal polynomials for the measure a 
giving rise to certain tridiagonal matrices (Jacobi matrices) so that the eigen- 
values of the n-th principal submatrix coincide with the nodes {aj^}"^^ i.e., 
with the zeros of the n-th orthogonal polynomial. Furthermore, the weights 
{Aj}"^^^ can be easily expressed in terms of the first component of the nor- 
malized eigenvectors. 

In this paper, we shall be concerned with the approximate calculation of inte- 
grals of 27r-periodic functions with respect to a positive measure /i on [— vr, vr] or 
more generally with integrals on the unit circle like I^{f) = J^^ f (e"^^^ dn{6). 
Here we will also propose as an estimation for //i(/) an n-point quadrature rule 
-^n(/) = jyj=i^jf{^j) with distinct nodes on the unit circle but now impos- 
ing exactness not for algebraic polynomials but trigonometric polynomials or 
more generally Laurent polynomials or functions of the form L{z) = J2'j=p C(jZ\ 
aj E C, p and q integers with p < q. Now, it should be recalled that Laurent 
polynomials on the real line were used by Jones and Thron in the early 1980 
in connection with continued fractions and strong moment problems (see [27] 
and [29]) and also implicitly in [30]. Their study, not only suffered a rapid 
development in the last decades giving rise to a theory of orthogonal Laurent 
polynomials on the real line (see e.g. [8], [14], [24], [28], [33] and [36]), but it 
was extended to an ampler context leading to a general theory of orthogonal 
rational functions (see [2]). 

On the other hand, the rapidly growing interest on problems on the unit circle, 
like quadratures, Szego polynomials and the trigonometric moment problem 
has suggested to develop a theory of orthogonal Laurent polynomials on the 
unit circle introduced by Thron in [36], continued in [26], [21], [11] and where 
the recent contributions of Cantero, Moral and Velazquez in [4] , [3] and [6] has 
meant an important and definitive impulse for the spectral analysis of certain 
problems on the unit circle. Here, it should be remarked that the theory of 
orthogonal Laurent polynomials on the unit circle establishes features totally 
different to the theory on the real line because of the close relation between 
orthogonal Laurent polynomials and the orthogonal polynomials on the unit 
circle (see [9]). 
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The purpose of this paper is to study orthogonal Laurent polynomials as well 
as the analysis and computation of the nodes and weights of the so-called 
Szego quadrature formulas and it is organized as follows: In Section 2, se- 
quences of orthogonal Laurent polynomials on the unit circle are constructed. 
They satisfy certain recurrence relations which are proved to be equivalent to 
the recurrences satisfied by the family of Szego polynomials, as shown in Sec- 
tion 3. The multiplication operator in the space of Laurent polynomials with 
a general ordering previously fixed is considered in Section 4. This operator 
plays on the unit circle the fundamental role in the five-diagonal represen- 
tation obtained in [3] analogous to the Jacobi matrices in the real line. Our 
main result of this section is to prove that this is the minimal representation 
from a different point of view than in [6]. In Section 5 a matrix approach 
to Szego quadrature formulas in the more natural framework of orthogonal 
Laurent polynomials on the unit circle is analyzed. Finally, we present some 
illustrative numerical examples on computation of the nodes and weights of 
these quadrature formulas, by way five-diagonal matrices versus Hcsscnberg 
matrices considering the so-called Rogers-Szego weight function in Section 6. 



2 Orthogonal Laurent polynomials on the unit circle. Preliminciry 
results 

We start this section with some convention for notations and some prelimi- 
nary results. We denote by T := {2; e C : \z\ = 1} and 13 := {z : \z\ < 1} 
the unit circle and the open disk on the complex plane respectively. P = C[z] 
is the complex vector space of polynomials in the variable z with complex 
coefficients, P„ := span{l, z,z^, . . . , 2;"} is the corresponding vector subspace 
of polynomials with degree less or equal than n while P_i = {0} is the triv- 
ial subspace. A := Cfz, 2;"^] denotes the complex vector subspace of Laurent 
polynomials in the variable z and for m, n G Z, m < n, we define the vector 
subspace Am.n = span{z'"\ z"^~^^ , . . . , z"-} . Also, for a given function f{z) we 
define the "substar-conjugate" as = f (l/z) whereas for a polynomial 

P{z) e P„ its reversed (or reciprocal) as P*{z) — z^P^{z) — z^P (1/^). 

Throughout the paper, we shall be dealing with a positive Borel measure /i 
supported on the unit circle T, normalized by the condition j^^ dfi{6) = 1 (i.e, 
a probability measure). As usual, the inner product induced by n{9) is given 

by 

if, 9),= j f {e^')JW)dm- 

— TT 

For our purposes, we start constructing a sequence of subspaces of Laurent 
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polynomials {>C„}^q satisfying 

dim {Cn) = n + 1 , £„ C £„+i , n = 0, 1, . . . . 

This can be done, by taking a sequence {p(^)}^o nonnegative integers 
such that p(0) = 0, < p(n) < n and s{n) = p{n) — p{n — 1) G {0, 1} for 
n = 1,2,.... In the sequel, a sequence {p(n)}^o satisfying these requirements 
will be said a "generating sequence" . Then, set 

Cn = A-p(„),g(n) = span {z^ : < j < , q{n) := n - (1) 



Observe that {g(n)}^Q is also a generating sequence and that A = U^o if 
and only if lim„^ooP(n) = lini„-^oo ?(n) = oo- Moreover, 



Cn © span{^«('*+^)} s(n + 1) = 
Cn ® span{z-P("+^)} i/ s{n + 1) = 1 



In any case, we will say that {p(n)}^Q has induced an "ordering" in A. Some- 
times we will need to define p(— 1) = and hence s(0) = 0. Now, by ap- 
plying the Gram-Schmidt orthogonalization procedure to £„, an orthogonal 
basis {%1)q{z), . . . ,%l)n{z)} can be obtained. If we repeat the process for each 
n = 1, 2, . . ., a sequence {V'nl^jl^o Laurent polynomials can be obtained, 
satisfying 



tpniz) e Cn\Cn-l , U = 1,2, . . . , Ipoiz) =C^0 



ii^niz) , Ipmiz)) n — l^n^n,m i '^n > , 5n,', 



if n ^ m 

1 if n — m 



(2) 



{ipniz)}'i^=o will be called a "sequence of orthogonal Laurent polynomials for 

the measure and the generating sequence {p(n)}^Q". It should be noted 
that the orders considered by Thron in [36] ("balanced" situation), expanding 
A in the ordered basis 



Ao,o 



A-1,0 , A_i^i , A_2,i , A_2,2 



A_ 



■3,2 



and 



Ao,o ; Ao,i , A_i^i , A_i^ 



A_ 



2,2 



A_ 



2,3 ; 



corresponds to p{n) = E 



n+l 



and p{n) = E ^ respectively, where as usual. 



E[x] denotes the integer part of x (see [3], [9] and [10] for other properties 
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for these particular orderings). In the sequel we will denote by {0„(z)}^q 
the sequence of monic orthogonal Laurent polynomials for the measure ^ and 
the generating sequence {p(n)}^0) that is, when the leading coefficients are 
equal to 1 for all n > (coefficients of z^^"^ or z"^*^") when s(n) = or 
= 1 respectively). Moreover, we will denote by {xn(-2)}^o sequence 
of orthonormal Laurent polynomials for the measure /x and the generating 
sequence {p(n)}^Q, i.e. when k,„ = 1 for all n > in (2). This sequence is 
also uniquely determined by assuming that the leading coefficient in Xni^) is 
positive for each n > 0. 

On the other hand, when taking p{n) = for all n = 0, 1, . . . then = 
^o,n = SO that, the n-th monic orthogonal Laurent polynomial coincides 
with the n-th monic Szego polynomial (see e.g. [35]) which will be denoted by 
Pn{z) for n = 0,1, .. .. This means means that po{z) = 1 and for each n > 1, 
Pn{z) e P„\P„_i is monic and satisfies 

{Pn{z), z')^ = (p;(z), z*)^ = 0, s = 0, 1, . . . , n - 1 , t^l,2,...,n 
{pn{z),z-),^{p:{z),l),>0. 



Moreover, we will denote by {v^ri(-2)}^o ^^'^ sequence of orthonormal polyno- 
mials on the unit circle for p{6), i.e., satisfying || ipn{z) ||^= {ipn{z), '^n{z))]l'^ ~ 
1 for all n > 0. This family is uniquely determined by assuming that the 
leading coefficient in </?„(z) is positive for each n > and they are related 
with the family of monic orthogonal polynomials by po{z) = (po{z) = 1 and 
Pn{z) = ln^n{z) with = {p^i^z) , Pn^z))]/"^ for aU n > 1. 

Explicit expressions for Szego polynomials are in general not available and 
in order to compute them we can make use of the following (Szego) forward 
recurrence relations (see e.g. [35]): 

Po{z) = p*o{z) = 1 

Pn{z) = ZPn-l{z) + SnPn-l{z) n>l (4) 
Pni^) = S^ZPn-l(z) + Pl_^{z) n > 1 



where So = 1 and 6n '■= Pn(0) for all n — 1,2,... are the so-called Schur 
parameters {Szego, reflection, Verhlunsky or Geronimus parameters, see [34]) 
with respect to p{6). Since the zeros of Pn{z) lie in D, they satisfy < 1 for 
n > 1. Now, if we introduce the sequence of nonnegative real numbers {r]n}'^=i 

by 

?7n = yr^^ e (0, 1] , n = l,2,..., (5) 
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then, a straightforward computation from (4) yields 




{pn-l{z),Pn-l{z)) 



and so, a forward recurrence for the family of orthonormal Szego polynomials 
is given by: 



nnVli^) = 5nZ^n-l{z) + ^l_i{z) 71 > 1. 

We conclude this section considering the following results proved in [3] and 
[11]. The first one establishes the relation between the families of orthogo- 
nal Laurent polynomials with respect to the generating sequences {p(n)}^Q 
and {q{n)}^^Q whereas the second one states the relation between the family 
of orthonormal and monic orthogonal Laurent polynomials for the generating 
sequence {p(n)}^o the family of orthonormal and Szego ordinary polyno- 
mials. This last one explains how to construct orthogonal Laurent polynomials 
on the unit circle from the sequence of Szego polynomials. Here it should be re- 
marked that the situation on the real line, i.e. when dealing with sequences of 
orthogonal Laurent polynomials with respect to a positive measure supported 
on the real line, is totally different (for details, see e.g. [14]). 

Proposition 2.1 Let {^n{z)}'^=Q be a sequence of orthogonal Laurent poly- 
nomials for the measure p and the generating sequence {q{n)}'^^Q. Then, 
in{z) = ipn*{z) for all n > 0, {'?/'n(^)}^o being a sequence of orthogonal Lau- 
rent polynomials for the measure p and the generating sequence {p(n)}^Q, 
where p{n) — n — q{n). 



Proposition 2.2 The families {0n(-2^)}^o {Xn{z)}'^^=o are the respective 
sequences of monic orthogonal and orthonormal Laurent polynomials on the 
unit circle for a measure p and the ordering induced by the generating sequence 



Vn^niz) = Zipn-l{z) + 5n(Pn-l{z) n>l 



(6) 



□ 



(Pn{z) 



and Xn{z) = < 



(7) 



□ 
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3 Recurrence relations 



We start with the next result which establishes a three-term recurrence relation 
for the monic orthogonal and orthonormal families of Laurent polynomials for 
the measure /i and certain (balanced) generating sequences, 

Proposition 3.1 Consider the familes {0„(2;)}^o ^'^^ {^n(-2)}J^o of monic 
orthogonal Laurent polynomials for the measure /i and the generating sequences 



p(n) = E 



n+l 
2 



and p(n) — E ^ respectively. Set 



Sn if n is even 
Sn if n is odd 



Then, 



(Pniz) = (An + (l)n-i{z) + r^2_^^(-i)>„_2(^) , n>2 



(f)Q{z) = 1 , = 5i + 



ln{z) = {An + An_,z(-'r^') 4^n-l{z) + vLl^^-'^^"'4>n-2{z) , n > 2 



(f)o{z) = 1 , (f)i{z) ^ Si + z. 



(1) 



(2) 



□ 



These recurrences were initially proved by Thron in [36] in the context of 
continued fractions. An alternative proof is given in [10] making use of (4) 
and Proposition 2.2. We will see now that the recurrences (4) and (1) (the 
same with (2)) are in fact equivalent. 

Theorem 3.2 The recurrences (4) and (1) are equivalent. 

Proof.- As above said, it remains to show that from the relations (1) we deduce 
the Szego recurrence (4). Indeed, from Proposition 2.2, (1) and by taking 
"super-star" conjugation it follows for A; > 1: 



P2k{z) = ((^2fe + ^2k-lZ) P2fe-l(^) + Vlk-l^ P2k-2{z) 



(3) 



= <^2feP2fe-l(^) + ^ \^2k-lP*2k-l{^) + Vlk-lZp2k-2{z) 
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P2fe-l(^) = (^2ifc-l^ + (^2fc-2)p2fe-2(^) +r?ifc_2P2fc-3(^) 



52k-\Zp2k-2{z) + 52k-2p2k-2{z) + %fe-2P2fe-3 



P2fc(^) = {^2kZ + (^2fe-l) P2fe-l(^) + ^ifc-lP2fe-2(^) 



= (^2ifc^P2fc-l(^) + 52ifc-lP2fe-l(2) +77ife_iP2fe-2(^) 



(4) 



(5) 



P2k-l{z) = (52ifc-l + 52k-2z) P*2k-2{^) + mk-2^ P2k~i{z) 



= 52fc-lP2fe-2(^) + ^ k2fe-2P2fc-2(^) + ^ife-2^P2fc-3 (^) 



(6) 



Clearly, the proof will be completed if from (3)-(6) we deduce that: 

P2k-l{z) = 62k-lP2k-li^) + Vlk-l^p2k-2{z) 



P2k-2{z) = <^2fc-2P2fc-2(^) +r72fc_2P2fc-3(^) 
P*2k-l{z) = (^2it-lP2fe-l(^) +?72fe-lP2fe-2(^) 
P2fe-2(^) = 52fc-2P2fe-2(^) +Vlk-2zp2k-^{z). 



(7) 
(8) 
(9) 

(10) 



For this purpose, it will be enough to check that (7) is valid since the proof of 
(8)- (10) follows in a similar way. Thus, set 

R2k-l{z) = 52k-lP*2k-l{^) + 'n2k-l^P2k-2{z) e P2fe-1. 

Hence, R2k-i{z) = OijPj{z). By comparision of monomials z^''~^ it fol- 

lows that a2k-i = 1- Since 

aj = {R2k-i{z),Pjiz))^ = S2k-i{P2k-ii^),Pjiz))n + vlk-i{zP2k-2iz),Pjiz))^, 

it follows from the orthogonality conditions for the Szego polynomials that 
aj = for all j = 1, ... ,2k — 2. Hence, R2k-i{z) = P2k-i{z) + cto- Finally, by 
comparision of monomials 2;° = 1 it follows that ao — 0. 

□ 



Moreover, the following recurrence relations were proved in [3] : 
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Proposition 3.3 The family of orthonormal Laurent polynomials for the mea- 
sure II and the generating sequence p{n) — E 



satisties 



^ X2n-i{z) 



X2n{z) 



-^2n-l(^2n —^2n-1^2n 



\ ^2n-l'72n 52n-1^72n / 



X2n-2{Z) 



\X2n-l\ 



+ 



-^2n(^2n+l ?72n?72n+l \ { X2n{z) 



^—52n^2n-\-l ^2nV2n+l J \X2n+l{z) 



, n>l 



□ 



A similar matrix recurrence is deduced if the generating sequence p{n) 



E 



n+l 
2 



is considered (see [3]). Once we have proved that the Szego recurrence 
(4) is equivalent with both given in Proposition 3.1 we will see now the equiv- 
alence with the recurrence given in Proposition 3.3. It will be proved in the 

case p(n) 



E 



ra+l 



IS 



. The proof when the generating sequence p{n) — E 
considered can "be established from Proposition 2.1 just by taking "super-star" 
conjugation. 

Theorem 3.4 The recurrence relations given in Propositions 3.1 and 3.3 are 
equivalent. 

Proof.- First of all, when the family of orthonormal Laurent polynomials is 
considered, (2) becomes 



Vn ■ Xn{z) = + An-lZ^ ^^"^') Xn-l{z) + 7]n-lZ^ ^^"^^ Xn-2{z) , 



n>2 



(12) 



Xi{^) 



i = 0,1. 



Since the recurrence (12) is equivalent to (4) and the recurrences given in 
Proposition 3.3 were obtained from this, it remains to prove that (11) implies 
(12). Prom Proposition 2.2 and since the generating sequence p(n) = ii^ | is 
considered, it follows that (11) is equivalent to the relations 

Z^ip2n-l{z) = ■n2nTl2n+lV2n+l{z) - V2j2n+l'P*2n(^) 

-S2n-lS2nZ^2n-l{z) - V2n-lS2nZ^2n-2{z) 
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^<^2n(^) = S2nV2n+1^2n+liz) - 52„52n+l¥'2n(^) ^^^^ 
+S2n-ir]2nZ^P2n-liz) + ?72n-l??2n^<^2n-2 (^) • 

Again, from Proposition 2.2 we have to prove from (13) and (14) that 

V2n^*2n{z) = {S2nZ + 52„-l) (p2n-l{z) + 7?2n-l'^2n-2 (^) (15) 

and 

ri2n+lV2n+l{z) = (52n+l + (^2n^) </'2n(^) + ?72n^^</?2n-l (^) (16) 

holds. By one hand, since r]2n 7^ it follows from (13) 

r]2nZ'^V2n-l{z) = (1 - |(^2nP) ^2n+l</^2n+l (^) - (1 - |(^2nP) (^2n+l<^2n(^) 
-?72n(^2n(^2n-l2;</'2n-l(^) " V2nV2n-lS2nZ'pln-2{z) 
= ^2n+l<^2n+l(^) " 52n+l¥'2n(^) " ^2n[52n^2n+l ¥5271+1 (^) 

-(^2nWl</'2n(^) + ?72n(^2n-l^</'2n-l (^) + V2nV2n-lZ<P*2n-2{^)] 

and now from (14) the relation (16) is deduced. By other hand, from (14) and 
since r]2n 7^ it follows that 

V2nZ(P2n{z) = S2nV2nV2n+l<P2n+l{z) - S2j2n+imnV*2n(^) + ^2n-lVL^V2n-l{z) 

+vLV2n-lZ(p*2ri-2iz) 

= (^2n?72n [?72n+l</^2n+l - 52n+l<^2n(^)] + 

S2n-l'P2n-l " (^2n-l |(^2n| V2n-1 + ?72n-l</'2n-2(^) " ^2n-l|(^2nP</?2n-2(^) ■ 
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Now, from (16) 



^2n<^2n(^) = ^2nr]2n [ri2nZ^2n-l{z) + 52n<^2n(^)] + <52n-l¥'2n-l " (52n-l | (52n P<^2n-1 
+V2n-l^*2n-2{z) " ^2n-l |(^2n| V2n-2(^) 

which is equivalent to 



V2n<P2ni^) " (^2n^</'2n-l (^) " (^2n-l</'2n-l(^) " ?72n-l</'2n-2 (^) 

So, (15) is proved since 772^ 7^ 0. 



□ 



Remark 3.5 Observe that the three-term recurrences given in Proposition 3.1 
involves multiplication by z and z^^ together whereas the recurrences given in 
Poposition 3.3 involves only multiplication by z. These last relations will play 
a fundamental role in the next section. 

Untill the present moment we have deduced recurrences for the families of or- 
thogonal Laurent polynomials when the generating sequences associated with 
the balanced ordcrings are considered. In the rest of this section we will con- 
sider an arbitrary generating sequence, starting with a three-term recurrence 
relation for {0„(2;)}^q involving multiplication by z and z~^ . Here, we define 
p(-l) = and so s(0) = 0. 

Theorem 3.6 The family of monic orthogonal Laurent polynomials {0„(2;)}^o 
with respect to the measure fi and the generating sequence {p(n)}^o satisfies 
for n >2 the three-term recurrence relation 



bn{z) = (AnBr, + CnZ'-''^'^^) 0„-l(^)+ ^^^^ 
(_l)l+.(n-2)-.(n-l)^^^^^2_^^l-.(„)-.(n-2)^^_^(^) 



with initial conditions 

0o(;^) = l , 0i(^) = Xi + z^-'^(^) 
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where Ki 



Si if s{l) = 
Titf s(l) = l 



and for n>2, 



if s{n) 7^ s{n — 1) 



c2s{n)-l 
"n-1 



if s{n) — s{n — 1) , s{n — 2) — , 



l-2s(n) 



if s(n) — s(n — 1) , s(n — 2) — 1 



(19) 



En 



5n if s{n) = 

) 

Sn if s{n) = 1 

1 if s{n) = s{n — 1) 

^Hn-D-sin) ^^^^ _^ _ , s(n - 2) = , 

_ (^)'^"^"'^""'^ if s{n) ^ s{n - 1) , s(n - 2) = 1 

5n if s{n — 1) = s{n) — 
5n if s{n — 1) = s{n) — 1 , 
1 if s{n — 1) 7^ s{n) 

l/5n-i if s{n - 2) = s(n - 1) = 



1/Sn-i if s{n - 2) = s(n - 1) = 1 . 
1 if s{n - 2) ^ s(n - 1) 



(20) 



(21) 



(22) 



(23) 



In the cases s{n — 2) = s(n — 1) the three-term recurrence relation holds if and 
only if5n-i ^ 0. 

Proof- The initial conditions (18) follows from Proposition 2.2 and (4). Con- 
sider, for n > 2, the eight cases (s(n), s(n — 1), s(n — 2)) = {i,j,k) with 
i, J, k e {0, 1}. As we said, the proofs in the balanced situations (0, 1, 0) and 
(1, 0, 1) were given in [10] making use of Proposition 2.2 and the Szego recur- 
rence (4). Proceeding in the same way in the remaining cases we obtain: 

(1) Case (0,0,0): 

M^) = + ^] 0n-l(^) - -^vl-lZ<l>n-2{^). (24) 
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(2) Case (0,0,1): 

0n(^) = + 0n-l(^) + 5n^^_l0n-2(^) • (25) 

(3) Case (0,1,1): 

(Pniz) = [Sn + =^ ) (pn-liz) - =^r)l_i(f)n-2{z) ■ (26) 
V (>n-lj (>n-l 

(4) Case (1,0,0): 

Mz) = (C+ J^] (Pn-l{z) - ■^vLl<Pn~2{z). (27) 

(5) Case (1,1,0): 

<Pn{z) = (^SnSn-l + 0n-l(2;) + 6nVl-l<Pn-2{z) ■ (28) 

(6) Case (1,1,1): 

Mz) =(= + -] Miiz) - j^vl-,M2{z). (29) 

Finally, all these three-term recurrences can be expressed as in (17) along with 
(19)-(23) in a whole pattern. 

□ 

Remcirk 3.7 When the family of orthonormal Laurent polynomials {Xn{z)}'^=Q 
is considered, the relation (17) becomes 

TJnXniz) = [AnBn + a.^^-'^^")) Xn-l{z)+ ^^^^ 



with initial conditions Xi{z) — n^'^^^^ji for i = 0, 1. 

In order to complete the equivalences between the recurrence relations it re- 
mains to show that the recurrence (17)-(23) implies the Szego recurrence (4). 
The proof in the balanced situations is given in Theorem 3.2 and the proof 
when a general generating sequence is considered can be done with similar 
arguments starting from the corresponding recurrence (24)- (29). In the next 
result we consider a general generating sequence and prove a similar result 
without using orthogonality conditions so that Proposition 2.2 can not be 
used. 
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Theorem 3.8 Let {p{n)}n>o be a generating sequence and define p{—l) — 0. 
Consider {s{n)}n>o defined by s{n) = p(n) — p(n — 1) and an arbitrary given 
sequence of complex numbers {(5n}n>o 5q = \ and ^ 1 for all n >1. 
Suppose that 5n-i 7^ if s{n — 2) = s{n — l). Let {4>n{z)}n>o be the sequence of 
Laurent polynomials defined by the recurrence relation (17) along with (19)- 
(23) and with initials conditions (18). Then, 



(1) Foralln > 0, 0„( z) is a monic Laurent polynomial with (f)n{z) G C^Cn—i- 

N„{z) 



(2) Write (j)n{z) = with N^iz) G Pn and set Fo{z) = No{z) = 1 and 



, NJz) if s(n) = 
Fn{z) = ^ ^ ^ ^ , n > 1. (31) 

N;{z) if s{n)^l 

Then, {Fn{z)}'^^i satisfies the recurrence relation 

F„{z) = zF„_,{z) + dr,F*_,{z) 

, n>l . (32) 

F:(z) ^ 5r,zF^^,(z) + F;_,{z) 

Proof- (1) Proceeding by induction it is proved for all n > 1 that (pni^) ^ 
Cr\Cn-i and that the coefficients of monomials z^^"^ and z"**^"^ are equal to 
1 when s{n) — and s{n) — 1 respectively, that the coefficient of monomial 

equal to 5n when s(n) = and that the coefficient of monomial 
equal to (5„ when s{n) = 1. The claim holds for n = 1 from the initial conditions 
(18). Hence, the proof can be achieved by induction from (17) by splitting into 
the eight cases (s(n), s{n — 1), s{n — 2)) = (i, j, k), with i, j, k e {0, 1}. 



(2) From (18) it clearly follows that A^o(^) = 1 and 



Si-^z if s(l) = 
S^z + lif s(l) = 1 



implying that Fi{z) = zFo{z) + 6iFq{z). Suppose from the hypothesis of 
induction that 

F^_^{z) = zF^_,{z) + Sr,-,F;_,{z) ^^^^ 
F*_,{z) = S;:r,zFr,.2{z) + F*_,{z) ■ 



Prom (17) it follows for n > 2 that 

^p(n) — [^nJ^n + <-^n^ J ^p(n)-s(n) + 

(_l)l+s{n-2)-s(n-l)r) Er? ^ _gl-g(n)-s(n-2) .^"-aW ^ 

I ^) ^p(n)-[s(n)+3(n-l)J 
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implying the following recurrence relation for the family of ordinary polyno- 
mials {Nn{z)}n>2- 

(_l)l+^(n-2)-.(n-l)^^^^^2_^^l+.(n-l)-.(n-2)^^_^(^)_ 



Wc consider now the four cases (s(n), s{n — 1), s(n — 2)) = (0, j, k), with j, k £ 
{0, 1}. Then, from (33) and relation (34) it follows: 



Case (0,0,0): 



NJz) 



zK 



^zN„_,{z) + ^^ 



Nn-l{z) - 7]l_^zNn-2{z) 

zN,,_2{z) + 5n-iK_^{z) - rjl_,zNr,.2{z) 

\Sn-l\'zNn-2{z)+5n-lN*_^{z) 



= zN,,_l{z) + Sn [5n-lzNn-2(z) + N*_2{z) 

= zNn-i{z) +6nN*_^{z). 



Case (0,0,1): 



Nn{z) = [Sn-lSn + z) Nn-l{z) + 8nril_^Nn-2{z) 



zNn-l{z) + 5n [5n-lNn-l{z) + ril_^Nn-2{z) 



ZN^_^{Z) + 5„ [5„_i^A^*_2(^) + Nr,-2{z) 

zN,,_,(z) + 5nN*_,(z). 



Case (0,1,0): 



Nn{z) = {Sn + Sn-iz) N„^i{z) + r]l_^Z^ Nn_2{z) 

= SnNn-l{z) + Z [Sn-lNn-l{z) + ril_-^zNn-2{z) 

= SnNn-l{z) + Z [\Sn-l\^zNn-2{z) + Sn-lK_2{z) + r]^„_iZNn-2{z) 
= dnN„-i{z) + Z [zNn-2{^) + dn-lN*_2{z) 

^5nN,,_,(z) + zN:_,(z). 
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Case (0,1,1): 

Nn{z) = + iVn-l(^) - ^y^_,zNr,.2{z) 

= 5„A^„_l(z) + Zj^ tzA^*_2(z) + Nn-2{z) - vLlNn-2 

"n — 1 L 

= SnNr,.i{z) + 2;=^ [5;:r72;iV;;_2(2;) + \Sn-l\^Nn-2{z) 

"n — 1 L 

= SnNn-l{z) + Z [zN*_2{z) + 5„_i7V„_2(^) 



Finally, the proof in the four remaining cases follows from the above proofs 
and from Proposition 2.1 just by taking super-star conjugation. 

□ 

We conclude this section by considering a Favard-type Theorem for the monic 

sequence of Laurent polynomials {0„(^)}^o a general generating sequence 
{p(n)}^Q. A proof in the ordinary polynomial situation (i.e. p{n) = for all n) 
is given in [26] and an alternative simpler approach in [15] (see also [1]). A proof 
based on the techniques introduced in the Chihara's book [7] in the balanced 
situations p{ri) — E 
presented here basec 



and p(n) — E ^ is given in [9] . Anyway, the proof 
upon the above equivalences between recurrences is very 
much simpler, and still valid when a general generating sequence is considered. 
For our purposes, we briefly recall the concept of orthogonality with respect 
to a Hermitian linear functional. Indeed, let {/^nlJ^-oo ^e a complex sequence 
satisfying = /JZ^ for all n = 0, 1, 2, . . . and denote by /i the linear functional 
defined on A by 



-j -I ^ 



— oo < p < q < +00. 



3=P 



In terms of // we define a bilinear functional < •, • >^ on A x A by 

< L, M >^:= 11 (l{z)M{1/z)) , L,MeA. 



(35) 



Now (see e.g. [26]), it will be said that the functional fi is quasi-deflnite if and 
only if the principal submatrices of the infinite Toeplitz moment matrix asso- 
ciated with {A*n}^oo nonsingular and positive-definite if the determinants 
of these matrices are positive. Quasi-definiteness is a necessary and sufficient 
condition for the existence of a family of orthogonal Laurent polynomials with 
respect to the linear functional (35) in the sense that there exists a sequence 
{Rn{z)}'I^^Q of Laurent polynomials satisfying Rn{z) e n = 1,2, .. . 

and < Rn{z),Rm{z) >^= k^Sn^rn ; 7^ 0. On the other hand, if the linear 
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functional n is positive-definite then, the associated linear functional (35) is 
an inner product on A x A and it holds < Rn{z),R„i{z) >^= knC^^.m with 
kn > 0. When = I for all n = 0, 1,2, . . . {R^{z)}^^o will be called "or- 
thonormal" . 

Theorem 3.9 (Favard) Let {p{n)}n>o be a generating sequence and define 
p(— 1) = and the sequence {s(n)}„>o by s{n) = p{n) — p{n — 1) G {0,1}. 

Consider an arbitrary given sequence of complex numbers {5„}n>o with 6o = 1, 
\6n\ 7^ 1 for all n > 1 and set {A„}„>i defined by \n = 1 — \6n\^ . Suppose also 
the restriction Sn-i ^ in the cases s{n — 2) = s{n — 1). Let {0„(-2)}n>o be the 
sequence of Laurent polynomials defined for n>2 by the recurrence relation 



with initial conditions (18) and where An, Bn, Cn, and E.^ are complex 
constants given by (19)-(23). Then, for a fixed /xq G IR\{0} there exists a 
unique quasi-definite linear functional fj, such that /x(l) = /iq and {4>n{z)}'^=o 
is the sequence of monic orthogonal Laurent polynomials with respect to /i and 
the ordering induced in A by the generating sequence {p{n)}n>o- Furthermore, 
if we take /Iq > then /i is positive definite if and only if \dn\ < 1 for all 



Proof- Prom Theorem 3.8 one sees that for n > 1, (f)n{z) G Cn\Cn-i so that 
one can write (t>n{z) — with Nn{z) e P„. If we define now the sequence 
of ordinary polynomials {pn{z)}n>Q by 



then, the relations (32)-(31) become (4). Thus, making use of the Favard's 
theorem in the ordinary polynomial situation (see e.g. [15] or [26]) we see that 
the family {pn(^)}n>o represents the family of monic orthogonal polynomials 
(Szego polynomials) with respect to a unique quasi-definite linear functional 
^l with //(I) = fiQ and positive definite if and only if //q > and \5n\ < 1 for all 
n = 1, 2, . . .. Finally, proof follows from Proposition 2.2 that still clearly holds 
when the measure /i is replaced for a quasi-definite (positive) linear functional. 




(36) 



+ {-iy+<-^)-s(-^)DnEr^Xn-lZ 



l_,(n)-s(n-2)0^_^^^-) 



n = l,2,.... 




□ 
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4 The multiplication operator in A 



Through this section it will play a fundamental role the multiplication operator 
defined on A, namely 



M: A — > A 
L(z) zL(z) 

As we have seen, if we consider the sequence of orthonormal Laurent polyno- 
mials with respect to the measure n and the generating sequence p{n) = 
for all n > then the n-th orthonormal Laurent polynomial coincides with 
the n-th orthonormal Szego polynomial, for all n > 0. Since the operator M 
lets P invariant, taking {ipn{z)}'^^Q as a basis of P then the following matrix 
representation of the restriction of M to P is obtained (see e.g. [19], [22] or 
[34]): 



n{d) 



W 

^1,0 ^1,1 ^1,2 
^2,0 ^2,1 ^2,2 ^2,3 



V 



The elements of this matrix with Hessenberg structure are given by 



hi 



-5j5i+i riLj+i Vkif i = 0, 1, . . . , i - 1 , 

-SiSi+i if j = i , 

Vi+i if i = ^ + 1 • 



(2) 



If we consider now the generating sequence p{n) = E ^ it follows from the 
recurrence given in Proposition 3.3 the following five-diagonal matrix (CMV 
representation, see [34]) for the multiplication operator M (see [3]) wich can 
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be also expressed as a product of two tri-diagonal ones: 



-5i rii 

-r)i52 -5i52 -r)253 772773 

^i?72 (^i?72 -S2S3 S2r]3 

-773(^4 -S3S4 -^74(^5 V^Vb 

7/37/4 53774 -54^5 ^4^75 



V : : : 

U 

-(^2 7/2 

772 5^ 



774 5l 



V 



-5i 771 

771 5^ 

-^3 ^73 

7/3 5^ 

-(^5 775 



(3) 



Moreover, it is easy to check that the matrix representation when the gener- 



ating sequence p{ri) — E 



n+l 



is considered is C(5)' 



The aim of this section is to analyze the structure of the matrix representa- 
tion for the multiplication operator M when an arbitrary generating sequence 
is considered. We start with the following result which is a generalization 
of Proposition 2.4 in [3] (here the generating sequences associated with the 
balanced orderings are previously fixed): 

Theorem 4.1 Let {Xn{z)}'^=Q he the sequence of orthonormal Laurent poly- 
nomials for the measure /j and the generating sequence {p{n)}'^^Q and suppose 
that hm„^oo9(^) = 00. Then, for each n > there exists k = k{n) > 1 and 
t = t{n) > 1 such that zxn{z) G span{xn-t{z) ■, ■ ■ -^Xn+kiz)}, i.e. 



ZXn{z) 



n+k 

an,sXs{z) 



an,s = {zXniz),Xsiz))f,. 



s=n—t 



Moreover, k = k{n) and t = t{n) are defined as follows: 

(1) k — 1 if s{n + 1) = and otherwise k>2 is defined satisfying 5(71 = 
• • ■ = s(72 + A; - 1) = 1, s{n + /c) = 0. 

(2) t = 1 if s{n — 1) = 1 and otherwise t >2 is defined satisfying s{n — 1) = 
■ ■ ■ ^ s{n + 1 - t) ^ 0, s{n - t) ^ 1. 
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Proof - On the one hand, since Xn{z) £ >C„ then 

with k = k{n) > 1. Observe that the existence of k is guaranteed from the 
condition hm„^oog(n) = oo. On the other hand, since Xn{z) -L then 

with i = t{n) > 1. Since 2;x„(2;) C Cn+k and 2;x„(2;) ± Cn-i-t the proof 
follows. 

□ 

From Theorem 4.1 we analyze now the matrix representation for the operator 
M with respect to a generating sequence with minimal number of diagonals. 
As we have seen, the balanced orderings gives rise to five-diagonal matrices 
and we must ask if there exists other possible generating sequences which gives 
rise to a matrix representation with less or equal number of diagonals than 
five. In this respect we start remarking that a five-diagonal representation is 
obtained when it happens 

1. k{n) = 1, t{n) < 3 for all n 

2. k{n) < 2, t{n) < 2 for all n (4) 

3. k{n) < 3, t{n) — 1 for all n. 

Hence, we have the following considerations: 

(1) Lemma 4.2 A five-diagonal representation is not obtained «/{s(n)}„>i 
contains three or more consecutive zeros or ones. 

Proof.- From Theorem 4.1 it follows that a block (s(n), s(n + 1), s{n + 
2), s(n + 3), s{n + 4)) = (1, 0, 0, 0, 1) implies (fc(n), k{n + 1), k{n + 2)) = 
(1, 1, 1) with k{n + 3) > 2 and {t{n + l),t{n + 2),t{n + ?>),t{n + 4)) = 
(1, 2, 3. 4) whereas a block (s(^), s{n + 1), s{n + 2), s(n + 3), s(n + 4)) = 
(0, 1, 1, 1, 0) implies {k{n), k{n + 1), k{n + 2), k{n + 3)) = (4, 3, 2, 1) and 
t(n+l) > 2 with (t(n + 2),t(n + 3),t(n + 4)) = (1, 1, 1). The proof follows 
since the condition (4) fulfill in both situations. Prom this, the proof when 
{s(n)}„>i contains more than three consecutive zeros or ones is trivial. 

□ 

As a consequence of Lemma 4.2, the number of consecutive zeros or ones 
in the sequence {s(n)}„>i is at most 2. 

(2) If the number of consecutive zeros or ones is just one for all n then it 
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corresponds to the generating sequences 



p{n) — E 



, p{n) = E 



n + 1 



and the five-diagonal matrix representation C{6) and C{6)^ are obtained 
respectively. 

Now, let us concentrate what happens if two consecutive zeros or ones appears 
in the sequence {s(n)}„>i. Indeed, a block of the form 

s{n + 1), s{n + 2), s{n + 3)) = (1, 0, 0, 1) 

implies t{n + 2) = 2, k{n + 2) > 2, t{n + 3) = 3, k{n + 3) > 1 and hence a 
five-diagonal representation is not obtained since condition (4) does not hold. 
Suppose now a block of the form 

s{n + 1), s{n + 2), s{n + 3)) = (0, 1, 1, 0). 

If s{n — 1) = then t{n -|- 1) = 3 and k{n -|- 1) = 2 whereas if s{n — 1) = 1 
then t{n) = 1, k{n) = 3 and t{n + 1) = k{n + 1) = 2, implying in both cases 
that the condition (4) is not satisfied. Observe that this argument is valid for 
all n > 2, but it really holds for n > 0. Indeed, if we consider the generating 
sequences p{n) = E for all n > 2 with p(0) = p(l) = or p{n) = e\^^ 
for all n > 2 with p(0) = and = 1 then it is easy to check that a non-five 
diagonal matrix representation is obtained. Summarizing, we can enunciate: 

Theorem 4.3 The matrix representation for the multiplication operator M 



or p(n) = E 



n+l 



is a five-diagonal matrix if and only if p(n) = E 
Moreover, this representation is the narrowest one in the sense that any ma 
trix representation for another different generating sequence gives rise to a 
I- diagonal matrix representation with I > 6. 



□ 



Remark 4.4 The proof of Theorem 4-3 was obtained using orthogonality con- 
ditions. This result has been also recently proved in [6] by using operator theory 
techniques. 

Example 4.5 Suppose that we expand A in the ordered basis 

{1, z, z-\ z\ z-\ z-\ z\ z\ z-\ z-\ z\ z\ z-\ z-\ . . .}. 

Here {s(n)}„>i = {0,1,0,1,1,0,0,1,1,0,0,1,1,...} and lim„^oo q'(^) = oo. 
Observe that zxo{z) G span{xQ{z) , Xi{z)} and zxi{z) E span{xo{z) , Xii^) , X2iz) , Xsiz)} . 
From the sequence {s{n)}n>i and Theorem 4-i we construct the sequences 
{t{n)}n>2 = {2, 1, 2, 1, 1, 2, 3, 1, 1, . . .} and {k{n)}n>2 = {1, 3, 2, 1, 1, 3, 2, 1,1,.. .}. 
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These sequences indicate us the number of nonzero elements in each file of the 
matrix. Now, the coefficients of the matrix can be obtained in terms of the 
orthonormal Szego polynomials by using the formula given in Theorem 4-i 
and Proposition 2.2, resulting expressions of the form {z''f{z),g{z))^ where 
I > and f{z),g{z) e {(fin{z)}n>o U {(Pni^)}n>o- Some of them are explic- 
itly calculated in terms of the sequence of Schur parameters (see [34, Ch. 1]). 
The remainder quantities can he obtained from these formulas and from the 
relations 

which holds since ^^i^) — ^"'Wni^) when z & T. The result is the following 
matrix representation with seven-diagonal structure: 



{z'+"'-^iPn{z),ipUz))a 



, t & Z , n,m > 0, 



-6i 


m 































-S1S2 



























Vim 
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-SgrigSio 




























V9 


-SgSio 

































-5io(5ii 



V • • • • • ■ ; ; ; ; ; ■) 

Finally wc conclude this section considering briefly the inverse multiplication 
operator defined on A, namely 

N : A — ^ A 
L{z) \L{z) 

The first result is an analogous to Theorem 4.1 which can be proved in the 
same way: 
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Theorem 4.6 Let {xn{z)}'^=o sequence of orthonormal Laurent poly- 

nomials for the measure fi and the generating sequence {p{n)}'^^Q and suppose 
that lim„^oop(n) = oo. Then, for each n > there exists k = k{n) > 1 and 
t^t{n) > 1 such that lxn{z) G span{Xn-i{z), ■ ■ ■,Xn+ki^)}> 

-j^ n+k 

-Xn{z) = J2 ^n,sXs{z) , dn,s = Oj^- (5) 
Z ~ 
s=n—t 



Moreover, k = k{n) and t = t{n) are defined as follows: 

(1) k — 1 if s{n+ 1) = 1 and otherwise k>2 is defined satisfying s{n-\- 1) = 
■■■ = s{n + k-l) = Q,s{n + k) = 1. 

(2) t = \ if s{n — 1) = and otherwise i > 2 is defined satisfying s{n — 1) = 
• • • = s(n + 1 - f) = 1, s(n - f) = 0. 

□ 

Now, proceeding as before, a similar result to Theorem 4.3 can be deduced for 
the inverse multiplication operator N. Furthermore, from (5) it follows that 
the matrix representation of the operator N when dealing with the balanced 

generating sequence p{n) — E ^ is C{S)* = C{S) . 



5 Quadrature formulEis on the unit circle 



Throughout this section we shall be concerned with the estimation of integrals 
on T of the form, 

TT 

W) = J f{e^')dm- (1) 



As usual, estimations of I^{f) may be produced when replacing in (1), f{z) by 
an appropriate approximating (interpolating) function L{z) so that Iij,{L) can 
be now easily computed. It seems reasonable to choose as an approximation to 
f{z) in (1) some appropriate Laurent polynomial, because of the density of A 
in C(t) = {/ : T — >• C , / continuous} with respect to the uniform norm (see 
e.g. [13] and [35]). So, the so-called "quadrature formulas on the unit circle" 
arise. Indeed, given the integral Ifiif), by an n-point quadrature formula on 
T we mean an expression like 

n 

W) = Yl ^jfi^j) ' , i¥'j , Zj er , j = i,...,n, (2) 
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where the nodes {zjYl^i and the coefficients or weights {)^j}^=i are chosen so 
that /„(/) exactly integrates /^(/) in subspaces of A with dimension as large 
as possible i.e., In{L) = I^{L) for any L e ^-p,q with p and q nonnegative 
integers depending on n with sum as great as possible. If we first try with 
subspaces of the form A_p^p, it can be easily checked that there can not exist 
an n-point quadrature formula In{f) as (2) to be exact in h.-n,n- Hence, it 
holds that p < n — 1 . In [1 1] the following "necessary condition" on the nodal 
polynomial is proved: 

Theorem 5.1 For n>l, let In{f) = E"=i ^jf{zj) with Zj E T, j ^ 1, . . . ,n 
be exact in A_(„_i)^„_i, and set Pn{z) — tVj=i{z — Zj). Then, 

Pn{z) = Cn [Pn{z) + T„Pn(^)] > K\ =1 , C'n = (l + Tn^n) (3) 

□ 

Moreover, in [26] it is proved the following converse result (sufficient conditions 
on the nodal polynomial): 

Theorem 5.2 Let Pn{z) be a polynomial of degree n given by (3) (up to a 
multiplicative factor). Then, 

(1) Pn{z) has exactly n distinct zeros zi, . . . ,Zn onT. 

(2) There exist positive real numbers Ai, . . . , A„ such that 

n 

In{f)^j:Xjf{zj)^lM) , V/ e A_(„_i),„_i. (4) 

□ 

The quadrature formula In{f) given by (4), and earlier introduced in [26] was 
called an "n-point Szego quadrature formula", representing the analogue on 
the unit circle of the Gaussian formulas for intervals of the real axis. However, 
it must be remarked two big differences in this respect: the nodes are not the 
zeros of the n-th orthogonal polynomial with respect to p and the n-point 
Szego formula is exact in A_(„_i)^„_i, whose dimension is 2n — 1 instead of 2n. 
Observe that since the nodes are the zeros of the n-th para-orthogonal poly- 
nomial characterized by (3), an one-parameter family of quadrature formulas 
exact in A_(„_i)^„_i arises. 

On the other hand, starting from a generating sequence {p{n)}'^^Q and because 
of the fact that >C„_i is a Chebyshev system on T of dimension n (since ^ T), 
for n distinct nodes Zi, . . . ,Zn on T, parameters Ai, . . . , A„ can be uniquely 
determined so that, /„(L) = I^{L) for all L G Cn-i- In order to recover Szego 
formulas in the natural framework of the orthogonal Laurent polynomials on 
the unit circle, and inspired by the ordinary polynomial situation, we will deal 



24 



with subspaces of A of the form Ln^r* = ^-\p{n)+q{r)],[q{n)+p{r)] with dimension 
n + r + 1 (observe that Cn~i C CnC^). Hence, from above, < r < n — 1 and 
we might analyze how large r = r(n) can be taken. The first step should be 
to consider r — n — 1, but a negative answer is proved in [11]: 

Theorem 5.3 There cannot exist ann-point quadrature formula like (2) with 
nodes on T which is exact in for any given arbitrary generating 

sequence {p{n)}n>o- 

□ 

The second step is to consider r = n — 2. For this purpose, we set A(n) ~ 
p{n) — p{n — 2) e {0, 1, 2}. The results obtained in [11] are summarized in: 

Theorem 5.4 Let {xniz)}^^^ be the sequence of orthonormal Laurent polyno- 
mials with respect to the measure /i and the ordering induced by the generating 
sequence {p{n)}'^^Q. Suppose that A(n) = p{n) — p{n — 2) = 1 and consider 

Rn{z, U) = Cn [rinXn{z) + T„Xn-l(^)] (5) 



being C„ 7^ and 



u-5nif s{n) = 1 
u — Sn if s{n) — 



(6) 



with u & T and {^njJ^Lo sequence of Schur parameters associated with /i. 
Then, 

(1) Rn{z,u) has exactly n distinct zeros on T. 

(2) If zi, . . . , Zn are the zeros of Rn{z, u), then there exist positive numbers 
Ai, . . . , A„ such that 

= E A,/(z,) = y /) , V/e«(„_2)*. (7) 

(3) There cannot exist an n-point quadrature formula with nodes on T to be 
exact in ^/ A(n) = or X{n) — 1. 

□ 

Thus, under the assumption that \{n) = p{n) — p{n — 2) = 1, we see that 
CnCi^n-2)* = ^-{n-i),n-i- Thcrcforc, the quadrature rule given by (7) coincides 
with an n-point Szego quadrature formula for ii{6) and, taking into account 
that the solutions of the finite difference equation A(n) = p{n) — p{n — 2) = 1 
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for n > 2 are given by 



tf p(0)=p(l) = 
i/p(0)=0, = l ' 

we see as the natural balanced orderings earlier introduced by Thron in [36] are 
again recovered. Furthermore, they are the only ones which produce quadra- 
ture formulas with nodes on T with a maximal domain of validity. On the 
other hand, as we have seen in the Section 4, these orderings correspond with 
the narrowest matricial representation of a sequence of orthonormal Laurent 
polynomials. 

In order to complete the construction of such quadrature formulas we give 
expressions for the weights also proved in [11]: 

Theorem 5.5 Let {Xn{z)}'^=o be the sequence of orthonormal Laurent poly- 
nomials with respect to the measure /i and the ordering induced by a generating 
sequence {p{n)}'^^Q. Then, the weights {Xj}^^i for the quadrature formula (7) 
are given for j = 1, . . . ,n by, either 




E 



n+l 



or 



23? 



+ (pin) - q{n)) \xn{zj 



(9) 



where the nodes {zj}'j^i are the zeros of Rn{z,u) given by (5), or equivalently 
the zeros of Pn{z) in (3). 

□ 

Once quadrature formulas on T with a maximal domain of validity have been 
constructed (nodes from (3) or (5) and weights by (8) or (9)) it seems that 
the following should be done is their effective computation. In this respect, it 
should be noticed that, untill the present moment and as far as we know, nu- 
merical experiments with Szego quadrature have mostly involved either mea- 
sures whose sequences of Szego polynomials are explicitly known or these 
polynomials have been computed by Levinson's algorithm (see e.g. [31], [12] 
and [21]). In this respect, the zeros of (3) or (5) can be found by using any 
standard root finding method available in the literature (as for a specific pro- 
cedure concerning rational modifications of the Lebesgue measure see also 
[20]). Hence, in the rest of the section we shall revise some strategies to ef- 
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fectivcly compute the nodes {-Zj}"=i and weights {Xj}]=i for an n-point Szego 
quadrature formula: 



4(/) = EV(^ 



Indeed, for the nodal polynomial Pn{z) — IYj=i {z — Zj) given by (3) we can 
write from (4) (see [4]): 

Pn{z) = Pn{z,u) = Zpn-l{z) + Upl_-^{z) , \u\ = 1. 

Now, it can be easily checked that {zipQ^z) , . . . , z(pn~2{z) , —uip^_-j^{z)} is an 
orthonormal basis of P„_i which must be related to {lpo{z), . . . ,ipn-i{z)} by 
an unitary matrix Un- Setting for n > 0, e„ = {pn{z) , pn{z)) ^ (recall that 
Cn = Ylk=i Vk for n > 1 and that ipn{z) = ^^-) it is straight forward to check 
that Un is the following irreducible Hessenberg matrix: 



d2,l d2,2 d2,3 



\ 



\ dn,l dn,2 dn,3 ' ' ' dn,n J 



(10) 



where 



di,j = < 



i/ i < n - 1 , j <i 
Vi ifi<n-l,j^i + l 



Hence, we can write 



zipo{z) 

Z(pn-2{z) 

\-uip*n_^{z) ) 



Ur.. 



^n-2{z) 
\Vn-l{z) I 



Observe that, if we introduce the notation 7i„(5o, • • • , ^n) to indicate the prin- 
cipal submatrix of 'H{5) (given by (1)) of order n with entries {5k}1=Q-i then 
Un — Ti-ni^o, ■ ■ ■ , Sn-i, u), that is, the principal submatrix of H(5) of order n 
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where the Schur parameter 5„ is replaced hy u E T. Now, 



\ 



ZCPn-2{z) 

[Z(pn-l{z) - Zipn-l{z) - Uipl_^{z) J 



Thus, we have the identity 

zVn{z) = UnVn{z) + bn{z) 



Vn-2{z) 





Pn{z,u) 



:i2) 



where 

Vniz) = {^Q{z),^l{z),---,^n-l{z))^ ] bn{z) = (o, ■ ■ ■ , 0, ^^^^^^] 

with |m| = 1. From (12) one sees that any zero ^ of Pniz^u) is an eigenvalue 
of Un with associated eigenvector Ki(^). So, let z^ be a zero of Pn{z^u) for 
J = 1, . . . , n and consider the corresponding normalized eigenvector of C/„: 



ELo \'^k{zj 



1/2- 



Now, taking into account that for any 2; e T, 199^(2;) | = 1x^(^)1 , from (8) it 
follows that 



(13) 



If we write Wj^z^ — {qo,j, ■ ■ ■ , Qn-i,j) and select out the first components of 



1 /2 

both sides of (13), we obtain q^j = \j ipo{zj). But, since we are dealing with 
a probability measure (/ diJ,{9) — 1) then (po{z) = 1 and hence, 

= qlj , i = 1, 
In short, the following theorem has been proved 

Theorem 5.6 Let /«(/) be the n-th Szego quadrature formula (7). Then, 

(1) The nodes {zj}'^^^ are the eigenvalues of Un — HniSo, ■ ■ ■ , Sn-i, u) given 
by (10)- (11), for all u e T. 

(2) The weights {Aj}^^^ are given by the first component of the normalized 
eigenvectors. 
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Finally, wc show an alternative approach to the computation of a Szcgo 
quadrature formula (7) by using truncations of the five-diagonal matrix C{d). 
In the next result we will use the matrix C„(5i, . . . , Sn-i, u), that is, the n-th 
principal submatrix of C{5) of order n where the Schur parameter Sn is replaced 
hy u E T. The first part has been already deduced in [5] by using operator 
theory techniques meanwhile our proof here presented is based on the recur- 
rence relations satisfied by the family of orthonormal Laurent polynomials. 
Without loss of generality, we can fix the ordering induced by p{n) = E ^ 
(recall that the matrix representation associated with the ordering induced by 
p{n) = E [^] is C{6)^). 

Theorem 5.7 Let In{f) be the n-th Szego quadrature formula (7). Then, 

(1) The nodes {zj}'^^]^ are the eigenvalues ofCn{Si, . . . , Sn-i, u), for allu e T. 

(2) The weights {\j}'f= I are given by the first component of the normalized 
eigenvectors. 

Proof- As we have seen, the nodes are the zeros of Rn{z) given by (5)-(6). 
Setting Xn{z) = (xo(^); ■ ■ ■ ;Xn-i(^))^ if follows from (11) that 



where 



being 



(0,...,0,A(z)f ifn^2k 



Tn{z) 



(0, . . . , 0, B{z), C{z)y if n^2k + l 



A{z) = r]2k-l{u - 52k)X2k-2{z) + 82k-\{u - S2k)X2k-i{z)- 

V2kS2k+lX2k{z) + V2kV2k+lX2k+l{z) 

B{z) = r]2k{u - 52k+i)X2kiz) + ri2kri2k+iX2k+i{z) 

C{z) = d2k{u - S2k+l)X2kiz) + S2kri2k+lX2k+liz). 

Now, on the one hand, in the odd case (n = 2A; + 1) it follows from (5)- (6) 
that ^ 

Tn{z) = (O, . . . , 0, r]n-lRn(z), 5n-lRn(z)^ 

implying that Tn{z) = (0, . . . , 0)^ if and only if Rn{z) = 0. On the other hand, 
from (12) and since u ET it follows in the even case {n — 2k) that 

A{z) = {U - 62k) Zr]2kX2k{z) - (?>2kZ + ^2k-\) X2k-\{z) 
+52jt-l(M - hk)X2k-\{z) - rj2kS2k+lX2k{z) + 
r]2k [(Wl + ^2kZ) X2k{z) + 'n2kZX2k-l{z)\ 
— UzRn{z) 
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and hence it holds again that Tn{z) = (0, . . . , 0)^ if and only if Rn{z) = 0. 
The proof of the second part follows from the same arguments as in Theorem 
5.6. 

□ 



6 Numerical examples 



In order to numerically illustrate the results given in the previous section 
we expose the computation of nodes and weights of the Szego quadrature 
formulas considering the absolutely continuous measure defined on [— tt, tt] by 

dii{e) = uj{e)de with 



27r 



2^ log ( J) 



exp 



J = -0O 



, ge(0,i) 



The corresponding monic orthogonal polynomials are the so-called Rogcrs- 
Szego g-polynomials. Throughout this section we also fix the ordering induced 
by the generating sequence p{n) = E . An explicit expression for such 
polynomials is given in [34, Ch. 1] and so from Proposition 2.2 the following ex- 
plicit expression for the corresponding monic orthogonal Laurent polynomials 
is deduced: 



2k 
j + k 



q 2 z-' 



1 

2k + I 
k-j 



if n — 2k 



g''^2^ ifn = 2k+l 



(2) 



where, as usual, the g-binomial coefficients 



are defined by. 



(n), = , (0)^ = 1 



n 
3 



(l-g")-(l-q"-^+l) 



(3) 
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Now, writting 02fc(^) = Ej=-fe%^^ and 02fe+i(^) = Ej=-(jt+i) then the 
coefficients {aj}^^_j^ and can be recursively computed by 

1 - y^+'+i 

Ofe = 1 , dj = -ttj+i^/q ^ _ , -k <j <k-l 



1 - q^~i 

= 1 , = -h,^ ^ _ , -(A; + 1) < i < A; - 1. 



The following two figures show the zeros of (t)io{z) and taking g = 

0.1, 0.25, 0.5, 0.75 and 0.9 (as usual, the circles below represent the unit circle). 
The data were computed by using the Jenkins and Traub root-finding method 
(see [25]), appropriate in this case since the polynomials has real coefficients, 
and they are similar to those given in [34, Ch. 8] where the distribution of 
zeros of Szego polynomials for this measure is also considered: 



Figure 1 {zeros of (pio^z)) 



g = 0.1 g = 0.25 g = 0.5 




Figure 2 [zeros of (j)ii{z)) 
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q = 0.75 



q = 0.9 






From Figures 1 and 2 one can observe that the zeros of (pioiz) are located on 
the circle : \z\ = g^^^j whereas the zeros of are located on the circle 

: \z\ = g~^'^^|, in accordance with Proposition 2.2 and the Mazel-Geronimo- 
Hayes theorem (see [32]). 

Now, we compute the nodes and weights of the quadrature formula (7) also 
for n = 10 and q = 0.1,0.25,0.5,0.75 and 0.9. The nodes are the zeros of 
Rio{z,u) given by (5) (we take u = 1). The corresponding sequence of Schur 
parameters are given by Sn — (— l)"?? for all n = 1, 2, . . . (see [34]). As for 
the weights, we can make use of (9) taking into account that 

Xn[Z) = 



;i-g)---(i-9") 

with (j)n{z) explicitly given by (2). The results are displayed in Tables 1-5. 
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Table 1 



Table 2 



q = 0.1 


n = 10 


Nodes 


Weights 


-0.940400 ± 0.34007Z 


0.0459602 


-0.531157 ± 0.847273i 


0.0669775 


0.0668824 ± 0.997761i 


0.100057 


0.624424 ± 0.781086i 


0.133157 


0.955949 ± 0.293533i 


0.153848 



q = 0.25 


n = 10 


Nodes 


Weights 


-0.922051 ± 0.387069Z 


0.0195773 


-0.473103 ± 0.881007Z 


0.0466391 


0.119954 ±0.992779ii 


0.0947585 


0.650270 ± 0.759703i 


0.150362 


0.959239 ± 0.282596Z 


0.188665 



Table 3 



Table 4 



q = 0.5 


n = 10 


Nodes 


Weights 


-0.842988 ± 0.537932? 


0.00312009 


-0.333209 ± 0.942853i 


0.0207928 


0.234605 ± 0.97209H 


0.0737936 


0.703537 ± 0.710659i 


0.163017 


0.965879 ± 0.258994i 


0.239274 



q = 0.75 


n = 10 


Nodes 


Weights 


-0.517559 ± 0.855648i 


0.000196919 


0.00961854 ± 0.999954Z 


0.00541542 


0.467501 ± 0.883993i 


0.0439839 


0.801825 ±0.597559i 


0.158275 


0.977622 ± 0.210369i 


0.292128 



33 



Table 5 



q = 0.9 


n = 10 


Nodes 


Weights 


0.112467 ± 0.993655i 


0.0000221869 


0.475746 ± 0.879582i 


0.000173911 


0.734593 ± 0.678508i 


0.0276214 


0.904709 ± 0.426031i 


0.146351 


0.989421 ± 0.145076i 


0.324221 



Remcirk 6.1 From the above tables one sees that the weights corresponding 
to a pair of complex conjugate nodes are equal. This directly follows from (8) 
since in this case the coefficients of the n-th orthonormal Laurent polynomial 
are real. 

In order to check the effectiveness of Theorems 5.6 and 5.7 we also com- 
pute for n = 10, g = 0.1,0.25,0.5,0.75,0.9 and u = 1 the eigenvalues and 
the first component of the normalized eigenvectors of Cio(5i, . . . , ^g, 1) and 
7^10 ((5o, • • • , <^9, 1), yielding the nodes and weights of the quadrature formula 
(7). The computations were made by using a standard eigenvalue-finding 
method and the results were exactly the same ones as those obtained be- 
fore. Moreover, the results of the computations for n — 11 are also showed in 
the following tables: 
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Table 6 



Table 7 



q = 0.1 


n = 11 


Nodes 


Weights 


-1 


0.03873768 


-0.814432 ± 0.580259Z 


0.04723999 


-0.355944 ± 0.934507i 


0.06942034 


0.199254 ±0.979948i 


0.09817396 


0.68359 lb 0.729866i 


0.12485883 


0.963208 ± 0.268755Z 


0.14093668 


Table 8 


q = 0.5 


n = 11 


Xodcs 


W'cigiils 


-1 


0.00105402 


-0.681725 ±0.731608i 


0.00606822 


-0.178603 ±0.98392H 


0.02834884 


0.335240 ± 0.942133i 


0.08185750 


0.745084 ± 0.666970Z 


0.16052782 


0.970795 ± 0.239909Z 


0.22291291 



q = 0.25 


n = 11 


Nodes 


Weights 


-1 


0.01403402 


-0.778821 ± 0.627246Z 


0.00025263 


-0.301023 ± 0.953617?; 


0.05384075 


0.243118 ± 0.969997i 


0.09739723 


0.70384 ± 0.710359Z 


0.14346753 


0.965731 ± 0.259545Z 


0.17367845 


Table 9 


q = 0.75 


n=ll 


Xodcs 


W(-igiils 


-1 


0.00000331 


-0.400355 ± 0.91636i 


0.00047627 


0.0972837 ± 0.995257i 


0.00836960 


0.517853 ± 0.85547i 


0.05348059 


0.821232 ± 0.570595i 


0.16811708 


0.979848 ± 0.199745Z 


0.29160532 



Table 10 



q = 0.9 


n = 11 


Nodes 


Weights 


-1 


0.00000000058 


0.149722 ±0.988728i 


0.00005234 


0.498881 ± 0.866671i 


0.00328985 


0.746631 ± 0.665238Z 


0.04447931 


0.909098 ± 0.416581Z 


0.21730036 


0.989911 ± 0.141688Z 


0.46353159 



Here, it should be remarked that when computations for higher values of n 
are required, special eigenvalue-finding methods should be considered because 
of the error propagation. In this respect, appropriate procedures for product 
of matrices (see [37]) could be considered due to the factorization of C{S) as 
product of two tri-diagonal matrices. Moreover, in [3] it is remarked that: 

(1) The computational cost by using techniques for banded matrices in com- 
parision with techniques for Hessenberg matrices is reduced. 

(2) In the C{5) matrix, the Schur parameters appear in only finitely many 
elements and hence, every modification of a finite number of Schur pa- 
rameters induces a finitely dimensional perturbation, something that is 
not true in the Hessenberg matrix. 

Throughout this section the computations were made using MATHEMATICA 
software. 
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